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Figure 1: Lego plot for double diff'ractive dissociation. 

1 Introduction 

A large rapidity gap (LRG) in an hadronic or DIS final state is experimentally defined as a large 
gap in the rj — (p lego plot where no hadrons are produced. A schematic description of a LRG 
di-jet diffractive final state is shown in Fig. 1. Historically, LRG events where singled outpP[2] 
as a signature for Higgs production due to a. W — W fusion sub process in hadron-hadron 
collisions at exceedingly high energies. At presently available energies, LRG di-jet central 
production is observed in the Tevatron[3|[?j[3| and HERA^. These measurements provide a 
unique opportunity to assess the asymptotic short distance behavior of the hard diffractive sub 
process responsible for the LRG di-jet in the final state and check our ability to adequately 
calculate it utilizing perturbative QCD (pQCD) methods. 

The production of two hard jets separated by a LRG can be explained by either: 

1) A fluctuation in the rapidity distribution of an inelastic event. The probability for such 
a fluctuation is proportional to e~^y^^, where A^ = Vi — y2 and L denotes the value of the 
correlation length. In a LRG event Ay ^ L and, thus, the probability for such a fluctuation is 
very small. Or: 

2) A QCD colorless t-channel exchange. We denote this exchange as a "hard Pomeron". Ac- 
cordingly, we shall be interested in Fg, the ratio between the LRG and the inclusive cross 
sections for di-jet production. It was noted by Bjorken0 that we have to distinguish between 
the calculated pQCD ratio Fs and the actual measured ratio fgap, where the proportionality 
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Figure 2: Pictorial definition of f gap- P and G represent color singlet and octet respectively. 



factor is called the " survival probability of a LRG" 
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There have been several suggestions how to calculate the survival probability. 

1) A simple procedure[7j has been suggested by Gotsman, Levin and Maor (GLM) to calculate 
{\ S \ ). The procedure is based on the GLM modeljS] which incorporates the effects of s- 
channel unitarity in high energy soft scattering. In the following I shall present this approach 
in some detail. 

2) An alternative procedure has been suggested by the Durham group |Hj. It will be discussed 
briefly in Sec. 6. 

3) Another alternative proposal ^Ij, associates the LRG phenomena with the color evaporation 
model (CEM). Like the previous listed options it explains the observed low values of the survival 
probability using the interplay between hard and soft physics. I shall discuss it very briefly in 
Sec. 6. 

The survival probability is the probability of a given LRG not to be filled by debris (partons 
and/or hadrons). These debris can originate from the soft rescattering of the spectator partons, 
and/or from the gluon radiation emitted by partons taking part in the hard interaction 



(I 5(S, Ay) 1^) = (I 5,,e^(Ay) 1^) ■ (I Sspecis) |'). 
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(I Siyrem \ ) cau bc Calculated in pQCD|llj. Its estimate 0, at presently available energies, is 
considerably smaller than ( | Sgpec I ) • 



The need for a reliable calculation of the survival probability is evident once we note that 
the values of fgap, measured at the Tevatronj3]|l] and HERA|n|, are very small and energy 
dependent: 

fgap{^/s = 630 GeV) = 1.60 ± 0.20% (DO), and 2.70 ± 0.90% (CDF), (3) 

fgapiVs = 1800 GeV) = 0.60 ± 0.20% (DO), and 1.13 ± 0.16% (CDF), (4) 

fgap{V^= 150 Gel/) = 7.00 ±2.00% (ZEUS). (5) 

These numbers are considerably smaller, and reduce with increasing energy, in contradiction to 
the pQCD expectations. These observations can be understood only if (| S* | ) is rather small 
and energy dependent. Indeed, a direct measurement of the energy dependence of fgap by DOjH] 
gives 

fgapiV7s= 630 GeV) 

fga,iV-s = 1800 GeV) -^-^il-^' W 

to be compared with a ratio of 2.7 ± 0.6 obtained from the individual numbers just quoted 
(Eq. (jni) and Eq. Q). As stated, the purpose of this short review is to present a relatively 
simple, and parameter free, procedure to calculate the survival probability and correlate its 
small value with the properties of high energy soft scattering. I shall show that the Tevatron 
data is perfectly consistent with the HERA data and that the calculations yield a natural 
interpretation of LRG di-jets observed in different channels. An important consequence of this 
analysis is that it provides significant support for the need to implement s-channel unitarity 
corrections to soft hadronic interactions in the HERA - Tevatron energy range. 

The organization of this paper is as follows: In Sec. 2 I briefly review the roll of s-channel 
unitarity in high energy soft scattering and its implications. This is followed in Sec. 3 by 
a suggested model by GLMlH^, compatible with unitarity constraints, which is based on the 
eikonal approximation. The calculation of the survival probabilities in the eikonal approach are 
presented and discussed in Sec. 4. The procedure presented incorporates the soft scattering 
basic measurements which replaces the need for free parameters. A generalization to a multi 
channel rescattering model is briefly reviewed in Sec. 5. A short discussion is presented in Sec. 
6. 



2 s-Channel Unitarity in High Energy Soft Scattering 

High energy soft scattering is commonly described by the Regge-pole model^2]- The theory was 
introduced more than 40 years ago and was soon after followed by a very rich phenomenology. 
The two key ingredients of this approach are the Pomeron whose t dependent trajectory is 
given by 

ap{t) = ap{0) + a'pt, (7) 



and a leading Regge trajectory 

anit) = ai?,(0) + a'j^t. (8) 

Donnachie and Landshoff (DL) have vigorously promoted [T^ an appealing and very simple 
Regge parametrization for total and forward differential elastic hadron-hadron cross sections in 
which 

and 



dt \ dt J 1^ 



e"^'\ (10) 
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jf-j and atot are related through the optical theorem. The forward elastic high energy 
exponential slope is given by 

B,i = Bo + 2a'pln {-) . (11) 



The data is excellently fitted with universal parameters: 

1) A supercritical Pomeron with an intercept ap(0) = 1 + e, where e = 0.08 accounts for the 
high energy growing total cross sections. Its fitted^j slope value is a'p = 0.25 GeV~^. 

2) The low energy data is controlled by the leading Regge trajectory whose universal parameters 
aredB] V = 0.45 and a'p = IGeV-^. 

DL offer a global fit to all available hadron-hadron and photon-hadron total and elastic cross 
section data. Note, though, that in reality only pp and 7p reactions have attained high enough 
energies in which the Pomeron parameters can be tested. We also note that the comprehensive 
DL analysis is limited to the properties of the elastic amplitude with no reference to diffractive 
scattering. 

Elastic and diffraction dissociation are similar processes which have predominantly forward 
imaginary amplitudes corresponding to the exchange of vacuum quantum numbers in the t- 
channel. As such, both channels are dominated in the high energy limit by a Pomeron exchange 
and are expected, in a simple Regge model, to exhibit rather similar dependences on the 
kinematic variable. Indeed, in the triple Regge limit applied to high energy single diffraction 
(SD) we getUg 

where M is the diffracted mass, and Gppp is the triple Pomeron vertex coupling. The virtue 
of this formalism is that it provides a strong correlation between the energy dependences of 
<^tot, (^eh Bel, and the energy and M^ dependences of agd- 

The simple model just presented is bound, eventually, to violate s-channel unitarity since 
a el and asd grow with energy as s^^, modulu logarithmic corrections, while atot grows only as s^. 
The detailed theoretical problems at stake are easily identified in an impact b-space formalism 
which is outlined below. 



The elastic scattering amplitude is normalized so that 

dOel 



dt 
and 



7T\Us,t)\\ (13) 



atot = Anlmfeiis^O). (14) 

The elastic amplitude in b-space is defined as 

aei{s,b) = ^Jdqe-'^-''Us,t), (15) 

where t = — g^. In this representation 

atot = 2 / dhlmaei{s,b), (16) 

aei = dh\ aei{s,b) |^ . (17) 

The common formulation of s-channel unitarity implies that | aei{s,b) \< 1, which is the 
black bound. We also have the analyticity/crossing bound. The position of this bound depends 
on the mass of the lightest exchange in the t-channel. Froissart|l(ij. at the time, considered 
it to be a vr meson. Our present thinking tends to consider it as the lightest glue-ball which 
can be exchanged in the t channel. The Froissart bound is obtained when the unitarity and 
analyticity/crossing bounds are combined. 

atot <Cln(^-^y (18) 

where C depends on the crossed channel exchanged mass. Note that both the unitarity and 
Froissart bounds are numerical, not functional. Consequently, the Froissart bound (with any 
reasonable value of C) is not applicable to present day accelerator physics, being much higher 
than the actual data. On the other hand, as we shall see, the unitarity bound is very relevant to 
the analysis of the presently available high energy data. A schematic illustration of the above 
is presented in Fig. 3. 

As noted, a simple Regge pole parametrization will eventually violate s-channel unitarity 
since the integrated elastic and diffractive cross sections grow with s much faster than the total 
cross section. The question is if this is a future problem to be confronted only at far higher 
energies than presently available, or is it a phenomena which can be identified through exper- 
imental signatures observed within the available high energy data base. It is an easy exercise 
to check that the DL modeljin], with its fitted global parameters, will violate the unitarity 
black limit for small b just above the present Tevatron energy. Indeed, CDF reports ^Zj that 
aei{b = 0, i/i = 1800) = 0.96 ±0.04. In a DL type model this should be refiected in e acquiring 
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Figure 3: A 'pictorial illustration of a high energy b-space elastic amplitude bounded by unitarity 
and analyticity/ crossing. In the illustration we have an input amplitude which violates the 
unitarity bound and an output amplitude obtained after a unitarization procedure. 



an energy dependence rather than the fixed fitted value. This is not observed, as yet, in the 
total and elastic cross section available data. We note, though, that the energy dependence of 
the experimental SD cross section j5j in the ISR-Tevatron energy range is much weaker than 
the power dependences observed for a^i- This observation is dramatically demonstrated once 
we check the ratio between f -^ j and f -^ j which goes down with energy. 



3 The GLM Screening Correction Model 

The theoretical difficulties, pointed out in the previous section, are eliminated once we take 
into account the corrections necessitated by unitarity. The problem is that enforcing unitarity 
is a model dependent procedure. In the following I shall confine myself to the GLM screening 
corrections (SC) modeljSj. In this semi realistic model the unitarity enforced SC are calculated 
utilizing a Glauber type eikonal approximation in which the scattering matrix is diagonal and 
only repeated elastic rescatterings are summed. Accordingly, we write 

a,i{s,h) = z (l - e-^(^'^)/2) , (19) 

where the opacity ^2 is a real function which is identified with the value of the imaginary part 
of a DL input super critical Pomeron. Analyticity and crossing symmetry are easily restored 
by the substitution s°^ -^ gap^-i-Kap/2 _ Combining Eq. ^, Eq. dTTj) and Eq. ^ we get 

a,ot = 2|d26(l-e~^(^'^)/2), (20) 

and 

a,„ = 2/rf2^(l-e-^(^'^)). (21) 



Since the scattering matrix is diagonal, the unitarity constraint is written as 

2Ima,i{s,h) = I a,i{sM |' + G^"(s,6). (22) 

Note that Eq. ^^ is a solution of Eq. ((221) • We also get 

G" = 1 - e-^^"'''). (23) 

The eikonal approximation can be summed analytically with a Gaussian input which corre- 
sponds to an exponential representation of the scattering amplitude in t-space (see Eq. (HU)))- 
This is, clearly, an over simplification, but it reproduces well the small t elastic cross section 
which contains the bulk of the data. 



n{s,h) = a{s^)Ts{s,h) = z/(s)e«^(^ (24) 



where 

and 

Rl(s) = ARl + Aa'pln{—). (26) 

So 

The soft profile is defined 



Some clarifications are in order so as to distinguish between the GLM model input and output: 

1) The power A, and accordingly z/, are input information and should not be confused with DL 
effective power e and its corresponding cross section. Obviously, A > e. 

2) The relation B^i = |-R|(s) depends mostly on the high b tail of the profile Ts- As such, SC 
do not change the output soft radius and this relation significantly. 

With this input we get 

atot = 2ttRI{s) f /n(^) + C - Eti-'^)] oc ln\s), (28) 



ain = vri?|(s) {ln[uis)] + C - Et[~uis)]) ex hn^is), (29) 

cTei = atot - cj,n = vri?|(s) (lui"^) +C - 2Ei{-^-^) + Ei[-v{s)]\ oc \ln\s), (30) 

where Ei{x) = J^^ ^dt, and C = 0.5773 is the Euler constant. 

The formalism just presented can be extended[Hj also to diffractive channels. The key 
observation is traced to Eq. (j^H|l from where we define 

Ps{s,b) = e-^(^''') (31) 

to be the probability that the two initial hadrons do not interact inelastically. i.e. there is no 
initial state rescattering in the inelastic sector at the given s, b values. Accordingly, the screened 
diffractive cross section is obtained by multiplying its b-space amplitude by this probability. In 
this approximation rescatterings in the diffractive final state are not included. Specifically, for 
SD we take the b-space transform of Eq. (fT^ multiplied by Ps{s, b) and obtain after integration 

r2/'c. ^ / o \ 2A 



-liM^ = 2^^m¥)\jp) ^""H^J ""^^'^^ [2.(s)^(--^^-) ^^""^''^^'''^^'^^' 

(32) 
where 

^^ {ip) = 2^° + ^° + '"'-^" {-i) ■ ('') 

8 



'"o ^ IGeV denotes the radius of the triple Pomeron vertex and can be safely neglected. Note 
1 

2^ 



that |i?| < i?| < i?|. In the ISR-Tevatron energy range i?| is just moderately larger than 



1 p2 
2^S- 



-(^.-^) - wi^fhwr '"' 



7(a, 2z/) denotes the incomplete Euler gamma function 



X 

a—1 „—z. 



7(a,x) = / z°--^e-^dz. (35) 

Jo 



In the high energy limit Eq. ()32|) simplifies to 



dM^ = 7r/gg(s)G'pppl — 1 oc/n(s). (36) 

The simple model, just presented, provides a significant improvement on the Regge DL 
model which serves as its input: 

1) It is compatible with unitarity (and Froissart) behavior in the high energy limit. In particular, 
the asymptotic ratio of ^ is kept at \ (Eq. (1^ - Ea. ^ '). 

2) Even though both atot and a^i change, at high enough energies, their dependence on s from a 
power to a ln'^{s) behavior, this change will be experimentally observed at an energy far higher 
than presently available. The model is, thus, compatible with DL in the ISR-Tevatron energy 
range. 

3) The SD cross section changes its power dependence on s to a ln{s) behavior. Checking the 
numerics of Eq. (j3(j|l this change is significant in the ISR-Tevatron energy range. Note that in 
the high energy limit ^^^^^ vanishes like 7-^. 

4) The significant difference between the energy dependence of the seemingly similar elastic 
and SD (and other diffractive) cross sections is understood once we examine their b-space 
amplitudes. The elastic amplitude is central (approximately given by a central Gaussian). Its 
DL parametrization implies that aei{s,b) will violate unitarity at small b somewhere between 
2 and 3 TeV cm. energy. Since this is confined to exceedingly small b values, it translates 
to very small changes in the normalization of the total and elastic cross sections. Screening is 
very different for the SD b-space amplitude. The screening suppression is obtained from the 
multiplication by Ps{s,b), of Eq. (jHT|) . which is very small at small b and approaches unity at 
higher b values. Whereas the screened elastic b-space amplitude maintains its centrality with 
a maximum at 6 = 0, the screened SD amplitude changes from a central to a peripheral b 
distribution, which is accompanied by a significant reduction of the output SD cross section. 
This phenomena is compatible with the Pumplin bound jTHj in which the sum of the elastic and 
diffractive b-space amplitudes is bounded by |. 

Regardless of its qualitative success, the toy GLM model does not provide a realistic repro- 
duction of the data. Its two important deficiencies are: 

9 



1) The model does not produce a good reproduction of the SD data [5]. This is amended in a 
multi channel eikonal model [7j in which diffractive intermediate rescatterings are also included. 

2) The model is applicable only for small t elastic scattering and it does not reproduce the high 
t diffractive dip experimentally observed at high energies. This is a consequence of our usage 
of a Gaussian b profile. To improve this deficiency we have to use a more elaborate profile [T^ . 



4 Survival Probability of LRG in the GLM Model 

Following Bjorken[2j, the survival probability, associated with the soft rescatterings of the 
spectator partons, is defined as the normalized integrated product of two quantities: 

, _ Jd%THis,b)Psis,b) 

1) The first element in Eq. ()37|1 is a convolution over the parton densities of the two interacting 
hadronic projectiles. This provides the uncorrected hard parton-parton collision cross section, 
leading, in our context, to di-jets. Following the discussion of soft scattering in Sec. 2, we 
assume a Gaussian hard profile 



H 



TiRii 



«!(-) 



(38) 



where Rh denotes the radius of the hard scattering process. This choice enables an analytical 
solution of Eq. ()37j] . More elaborate choices require a numerical evaluation of this equation. 
2) The second element is the probability Ps{s,b), defined in Eq. (jSH), that no inelastic soft 
interaction takes place between the incoming projectiles at impact parameter b and cm. energy 
square s. We define 

, , -R|(s) interaction area for soft collisions , . 

R'h{s) interaction area for hard collisions 

affi^s) grows logarithmically with s. As stated, Eq. (|H7j) can be analytically evaluated with our 
choice of Gaussian profiles and we get 

/, Q |2v _ aHis)j[aHis),u{s)] 

z/(s) was defined in Eq. (^3]) and the 7 function in Eq. ()35|) . 

A straight forward estimate of the LRG survival probabilities can be derived |2] [7j |9] [Toj 
from various models assumed to describe the soft and hard interactions. An alternative method 
suggested by GLMp| is based on the GLM eikonal model reviewed in Sec. 3. This calculation 
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of (I Sspec I ) is executed with no further recourse to the theoretical models, utilizing just the 
experimental data from which we can directly deduce the input needed to solve Eq. (j40|) . To 
this end we have to assess the values of Rj^, i?| and i^{s). 

The value of R% is estimated from two independent experiments: 

1) In analogy to the relation i?| ~ 2i?ei, we can estimate the value of the hard radius from the 
exponential slope of a clean hard process R^ ~ 2Bhard- We use the excellent HERA data [201 
on 7 + p — i> J/\l/ + p, where -Bj/* — 4:GeV^'^ with a small logarithmic dependence on the 
incoming energy. 

2) The above is compatible with the CDF estimate of the double parton cross section. This 
cross section is connected through factorization to the effective hard cross section 



(^DP 



m- 



<yincll{2jets) (Tincl2{2jets) 



^^eff 



(41) 



where the factor m is equal to 1 for identical pairs and 2 for different pairs, a^fj relates directly 
to the hard profile (Eq. (I^Hj)) 



1 



a 



H 

eff 



2t^R]j 



(fbTlib). 



(42) 



The experimental value [2Ij cr^ 



14.5 ± 1.7 ± 2.3 mb suggests that Rj^ = 5 - 7GeV-'^. 

8GeV~^. This is a conservative 
choice which may be slightly changed with the improvement of the Tevatron estimates. 



// 
In our calculations we took the HERA deduced value Rj^ 



The values of R's{s) are taken from the experimental elastic slope of the high energy pp 
data. The values of i/(s) can be obtained in the GLM model directly from the measured values 
of ^. Checking Eq. (J2HI) and Eq. (fHUj). We note that the (s dependent) ratio of elastic to total 
high energy cross sections provides explicit information on z^(s). Once we have determined z/(s) 
and anis), the survival probability is obtained from Eq. (j40p . 

The important result of our calculation, which is displayed in Fig. 4, is that (| Sspec I ) 
decreases as we go to higher energies. This is traced to the experimentally observed increase 
with energy of z/(s) and anis). To further illustrate this feature, we show a contour plot of 
(I Sspec I ) against u^s) and an^s) in Fig. 5. 

Our results for the survival probability at the reported energies: 



v^ (GeV) 


\ ^spec \ 1 


150 
630 
1800 


16.3% 
12.8% 
5.9% 



provide a qualitative explanation for the smallness and energy dependence observed experimen- 
tally. However, for a quantitative comparison between the reported fgap data (Eq. (jH)) - Eq. (jHI)) 
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Figure 4: (| Sgpec \ ) as a function of u{s) and anis) 
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Figure 5: A contour plot of (| Sgpec \ ) against uls) and a_f/(s) 
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and our results we have to incorporate a pQCD estimate of Fg. I shall return to this issue 
in the discussion. To bypass this difficulty it is constructive to check the ratio of the survival 
probabilities calculated at 630 and 1800 GeV 

^^^''-[^^^ = ^^^ =2.2 ±0.2. (43) 

\| ^spec I /^=1800 

Our calculated ratio is fully substantiated by the DO result j2] in which the corresponding 
survival probability (rather than the fgap) ratio was estimated at 2.2 ± 0.8. Note that for our 
calculations at 630 GeV we had to interpolate the input data information from 546 GeV. 



5 Extension to a Multi Channel Model 

In as much as the GLM results on the survival probability are very satisfactory, the over all 
quality of the GLM predictions is not that good. The GLM model |8, was originally conceived 
so as to explain the exceptional mild energy dependence of soft diffractive cross sections. It 
suggested that a high energy change from a power to a logarithmic energy dependence of the 
diffractive cross sections is a signature indicating the onsetting of s-channel unitarity corrections. 
The model actual calculations do, indeed, result in a qualitative significant difference between 
the energy dependence of the total and elastic cross sections, on the one hand, and diffractive 
cross sections, on the other hand. However, they do not provide an adequate reproduction of the 
SD cross sections|n] in the ISR-Tevatron energy range. A possible remedy to this deficiency of 
the GLM model is to replace the one channel eikonal model, in which only elastic rescatterings 
are included, with a multi channel eikonal model in which inelastic diffractive intermediate 
rescatterings are included as well. This is not a difficult technical problem j^l], but, in our 
context, we have to secure that a multi channel model does improve the diffractive (specifically 
SD) predictions of the GLM model, while maintaining, simultaneously, its excellent results on 
LRGtZj. 

The implicit assumptions of the simple GLM model are: 

1) Hadrons are the correct degrees of freedom at high energies. 

2) At high energy Rea^i ^ Iniaei- 

3) (Tdiff < (Tel- 

4) Only the fastest partons can interact with each other. 

Clearly, the last two assumptions are the weakest of the four. There is no question that at 
small enough x (high enough energy) one has to take into account (| Sbrem{s) \ ). However, an 
explicit pQCD calculation ^Tj has shown it to be rather small in the energy range of interest, 
i.e. ISR-Tevatron. Here, I present a re-formulation of the GLM model. The goal is to construct 
a multi channel eikonal model in which the rescattering can be either elastic or diffractive. 
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In the simplest approximation we consider diffraction as a single hadronic state. We have, 
thus, two orthogonal wave functions 

i^Hh I $,) = 0. (44) 

"^h is the wave function of the incoming hadron, and "$£> is the wave function of the outgoing 
diffractively produced particles, initiated by the incoming hadron. Denote the interaction 
operator by T and consider two wave functions \l/i and "$2 which are diagonal with respect to 
T. The amplitude of the interaction is given by 

A,fc = (*i^fc I T I "^i'^k') = 4> ki' ^k,k'- (45) 

In a 2 X 2 model i,k = 1, 2. The amplitude Afj^ satisfies the diagonal unitarity condition (see 
Eq. dUD) 

2JmAi(s,6) =Mi(.,6) p +G-,(s,6), (46) 

for which we write the solution 

Aj(.,6) =^[l-e--^y (47) 

and 

Qin^ = i_e-f^'.^(^'*). (48) 

^i,k{s, b) is the opacity of the {i, k) —th channel with a wave function \I'j x \I'jt. The probability 
that two hadronic states i and k do not interact inelastically is 

P.4s,b) = e-^^-'^^^'^) = {1 - A^Us,b)Y. (49) 

In this representation \E'/j and \E'rf can be written as 

<i/h = «^i + /5^2, (50) 

^d = -/5^i + a*2- (51) 

Since | \I//i ^ = 1, we have 

a'^ + p'^ = 1. (52) 

The wave function of the final state is 

^/ = I T I ^/, X ^,,) = a^Ai,i^i X ^1 + a/5Ai,2{^i x ^2 + ^2 x ^1} + /3^A2,2*2 x ^2- (53) 

Since Afj^ is a 2 x 2 matrix, we have to consider 4 possible rescattering processes. However, in 
the case of a pp (or pp) collision, single diffraction at the proton vertex equals single diffraction 
at the antiproton vertex, i.e. Ai 2 = v42,i and we end with 3 channels whose b-space amplitudes 
are given by 

aei{s,b) = {^h X ^h I */) = aMi,i + 2a2/5Mi,2 + /3M2,2, (54) 
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asd{s,b) = {^h X ^d I ^/) = «/5{«Mi,i + {a' - /?^)Ai,2 + /?M2,2}, (55) 

addis,b) = (^rf X ^rf I ^;) = a2;52|^^^^ _ ^a,^^ + ^2,2}- (56) 

Afi^ and G*'^ are given by Eq. (P7|) and Eq. (|iS|) . The input values of fij^ satisfy the Regge 
factorization property 

= "■'-';■'''"-;_''-■''> . (57) 

i2i,i(s,6) 

As in the single channel GLM model we simplify the calculation assuming a Gaussian b- 
space distribution of the input opacities Qi^k- 

^1,1 = i^iX, (58) 

^2,2 = iy2X'/^'-^\ (59) 

f^i,2 = ^'yM2 - r)X\ (60) 
X = e~^ /^i,i('') and r = R\ ^{s) / R\ ^{s) ■ The factorizable radii are given by 

Rlk = 2^0 + 2i?lo + 4a'p/n(-). (61) 

■50 

The opacity expressions just presented allow us to express all physical observable of interest 
as functions of i^i, z/2, r and (3. The first three variables depend on s while /3 is a constant of 
the model. The determination of these variables enables us to produce a global fit [22] to the 
total cross sections as well as the elastic, single and double diffractive integrated cross sections. 
This has been done in a two channel model, in which add is neglected and a more detailed 
three channel model. The main conclusion of these studies is that the extension of the GLM 
model to a multi channel eikonal results with a good overall reproduction of the data. The 
results maintain the b-space peripherality of the diffractive output amplitudes and satisfies 
the Pumplin bound. Note that since different experimental groups have been using different 
algorithms to define diffraction, the SD experimental points are too scattered to enable a tight 
theoretical reproduction of the data. For details, see Refs. j3], [S] and j22j and Fig. 6. 

To complete this discussion on the generalization of the GLM model we have to show that 
the multi channel model reproduces the excellent results we have obtained for the survival 
probability in the simple, single channel, model. This is not self evident. In the single channel 
model we have correlated the decrease with energy of (| Sgpec P) mainly with the corresponding 
power like increase of Rei = Cei/atot, see Fig. 7. In a 3 channel model Rei is replaced by Rd = 
((Tez + asd + CTddl/o'tot, which is approximately constant, see Fig. 8. Indeed, a calculation|^ 
suggests that 

{\Sspec?) = (1 - Rd)\ (62) 
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Figure 6: Integrated SD data and a 3 channel model fit^ 
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which yields a reasonable value of the survival probability, but fails to reproduce its energy 
dependence. Upon investigation, one finds that this problem originates from the neglect of the 
differential b-dependence of the various scattering amplitudes in Ref . J2S1 • From the 3 channel 
GLM model we obtain 

(I Sspec 1^) = -j^, (63) 



where 

N{s) = / —{(1 - /52)2Pi,iZ/iai,iX'^M + 2(1 - /?')/3'Pi,2^i,2ai,2X'^^'^'- (64) 

Jo X 



/o 'Y 



and 



"^dX 



Dis) = / —{(1 - /32)2,.,ai,iX«M + 2(1 - /3')/3'u,,,a,,,X^^'-^ (65) 

Jo X 

+ f3^U2a2,2X^'-''l{2 -r)}. 
We denote i/i_2 = -i/i^iZ/2(2 — r) and aj^^ = -^. 



A b-dependant calculation [7j [221 ^^ (I ^spec P) and its energy dependence reproduces the 
DO, CDF and HERA results [3] j3][n], with (3 ~ 0.5, which is a very reasonable value. 



6 Discussion 

Both the single channel and the improved multichannel GLM models, with which we calculate 
the survival probability and its energy dependence, do not attempt to calculate the pQCD 
ratio Fg which is an external input. Accordingly, the GLM model, on its own, cannot provide 
a calculation of fgap- Never the less, the model provides an important, and relatively simple, 
direct correlation between the soft scattering experimental features and (| Sgpec I )• i-e. given 
the soft scattering experimental features and the experimental value oi Bhard, we can predict the 
survival probabilities. As such even the partial experimental information we have on (| Sgpec \ ) 
from LRG di-jet production serves as an excellent probe of the roll of s-channel unitarity in 
high energy soft scattering, both elastic and diffractive. 

The GLM model is conceptually different from both the Durham|H] and the CEMp^ al- 
ternative approaches. All three models calculate (| Sspec I ) using the interplay between hard 
and soft physics and obtain similar values for the survival probabilities. However, we note 
that the GLM model does not dependent on input free parameters, at the cost of having a 
more limited predictive power than the other two models. At the present stage this is not a 
significant disadvantage since the experimental data base is rather small. Both Durham and 
GEM are partonic models which have a wider predictive power at the cost of depending on 
a rather extensive input. Note that even though both GLM and Durham are multichannel 
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models, they are dynamically different. GLM multi channel calculation relates to the diversity 
of the intermediate rescatterings, i.e. elastic and diffractive. In the Durham model the term 
"multi channel" relates to the input information. The model investigates two options: i) Small 
and large dipoles, ii) Valence and sea partons. 
Both options give similar results which are compatible with GLM. 

Much attention has been given recently to the compatibility between the Tevatron and 
DESY data. Clearly, rather than depending on a pdf input to calculate Fg we may use the 
gluon structure function inferred from the diffractive HERA datajH] and use it as input to the 
calculation of Fg at the Tevatron. This is an over simplified procedure ignoring the roll of the 
survival probability, and it has led to speculations about a possible breaking of QCD or Regge 
factorization or both. Once the dependence of (| Sgpec | ) on s and Q"^ is investigated it is an 
easy exercise to re establish the compatibility between the Tevatron data and DESY photo 
and DIS diffractive data. For a detailed discussion see Refs. jS] and p. In this context we 
note that the recent CDF data on multi LRG gap jet production[5] provides additional support 
to the philosophy advanced in this summary. As long as Rfj {multi LRG) ~ Rh{LRG) we 
anticipate the two survival probabilities to be approximately the same as, indeed, is reported 
experimentally. 
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